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Abstract  Powerful  nonlinear  approaches  for  missile  autopilot  design  have  recently 
emerged  in  the  literature,  which  have  the  potential  to  deliver  improved  missile  per¬ 
formance.  However,  the  lack  of  computational  methods  has  made  it  difficult  for  the 
practicing  engineers  to  exploit  these  techniques  in  routine  applications.  Another 
factor  that  has  slowed  their  application  is  that  the  missile  models  are  generally 
available  in  the  form  of  simulations,  rather  than  as  compact  set  of  differential- 
algebraic  equations.  This  paper  discusses  five  different  approaches  for  computer- 
aided  nonlinear  control  system  design  that  ameliorate  these  difficulties.  Since  these 
design  techniques  are  based  on  simulation  models,  they  enable  direct  synthesis  of 
nonlinear  autopilots  using  missile  models  of  arbitrary  complexity.  Airframe  stabi¬ 
lization  of  a  nonlinear,  longitudinal  missile  model  is  used  to  illustrate  the  design 
techniques. 


1  Introduction 

Methods  for  nonlinear  control  system  design  have  been  of  significant  interest  in 
the  recent  literature  [1  -  19].  By  enabling  the  design  of  missile  autopilots  with¬ 
out  employing  Taylor  series  linearization  and  subsequent  gain  scheduling,  these 
methods  have  the  potential  to  enhance  the  missile  performance.  While  some  of 
these  techniques  have  advanced  to  a  point  where  they  can  be  routinely  employed, 
the  nonlinear  design  processes  are  largely  based  on  algebraic  manipulations  of  the 
underlying  mathematical  model  of  the  system  to  be  controlled.  Although  con¬ 
troller  design  using  algebraic  manipulations  are  effective  in  simpler  problems,  it 
becomes  increasingly  onerous  to  employ  them  in  practical  situations  where  the 
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missile  model  may  contain  complex  nonlinearities  that  may  not  be  describable  in 
terms  of  symbolic  expressions,  such  as  sensor-actuator  nonlinearities  and  lookup 
tables. 

In  order  to  motivate  subsequent  discussions,  the  flight  control  system  of  a 
typical  homing  missile  is  illustrated  in  Figure  1.  The  state  estimator,  the  guidance 
law  and  the  autopilot  form  the  three  major  components  of  the  system.  The  state 
estimator  uses  measurements  from  onboard  sensors  together  with  a  set  of  missile- 
target  mathematical  models  to  generate  estimates  of  the  missile  and  the  target 
states.  The  guidance  law  uses  the  missile  and  target  state  estimates  to  generate 
commands  for  the  autopilot.  The  guidance  commands  are  typically  in  the  form  of 
lateral  acceleration  components  ay,azthat  will  achieve  target  interception. 

The  autopilot  has  the  responsibility  for  generating  actuator  inputs  for  tracking 
the  guidance  commands  while  stabilizing  the  missile  airframe.  In  addition  to  the 
identification  of  the  three  flight  control  subsystems,  Figure  1  also  illustrates  the  sets 
of  variables  employed  by  each  of  them.  For  instance,  the  state  estimator  uses  the 
missile  kinematics,  an  assumed  target  maneuver  model  and  measurements  from  the 
seeker  to  estimate  the  target  states.  The  seeker  measurements  may  consist  of  line 
of  sight  rates  cj,  line  of  sight  angles  a,  range  i?and  range  rate  R.  Depending  upon 
the  assumed  model  of  the  target,  the  state  estimates  may  include  target  relative 
position,  velocity  and  acceleration  vectors.  The  guidance  law  uses  estimated  target 
states  and  the  missile  position  and  velocity  vectors  to  generate  lateral  acceleration 
commands  for  the  autopilot.  The  autopilot  tracks  these  commands  while  ensuring 
the  stability  of  the  missile  short  period  dynamics  [20]  consisting  of  the  angle  of 
attack  a,  angle  of  sideslip  /?,  pitch  rate  gand  yaw  body  rate  rstates.  The  roll  rate 
pis  often  regulated  about  zero. 

Traditional  approach  to  autopilot  design  is  to  first  linearize  the  missile  short 
period  dynamics  about  an  operating  condition,  and  then  to  apply  linear  control 
theory  to  synthesize  a  feedback  controller.  This  process  is  repeated  at  multiple 
operating  conditions  and  the  controllers  are  then  scheduled  with  respect  to  the 
flight  conditions.  Typical  scheduling  variables  include  Mach  number,  dynamic 
pressure,  angle  of  attack  and  angle  of  sideslip.  In  some  missiles  capable  of  operating 
over  a  wide  range  of  altitudes,  the  controllers  may  also  be  scheduled  with  respect 
to  altitude.  While  linear  control  theory  provides  elegant  algorithms  for  controller 
synthesis,  the  gain  scheduling  step  is  largely  a  trial  and  error  process.  Often, 
the  gain  scheduling  process  can  consume  a  significant  portion  of  the  design  effort. 
Moreover,  the  stability  and  performance  guarantees  provided  by  the  linear  design 
techniques  are  generally  diluted  by  the  gain  scheduling  process. 

Nonlinear  control  system  design  methods  seek  to  eliminate  the  gain  schedul¬ 
ing  process,  without  compromising  the  performance  and  stability  properties  of 
the  closed-loop  system.  Moreover,  these  techniques  provide  a  natural  setting  for 
including  nonlinearities  such  as  actuator  saturation  in  the  design  process.  The  ob¬ 
jective  of  the  present  paper  is  to  advance  a  set  of  nonlinear  control  system  design 
techniques  that  allow  the  direct  use  of  numerical  simulation  models  of  the  missile 
for  nonlinear  autopilot  synthesis.  The  analyst  exercises  control  over  the  design 
processes  through  the  selection  of  structural  properties  of  the  controller  and  the 
parameters  that  govern  a  specific  design  technique.  Design  techniques  presented 
here  can  handle  a  large  class  of  system  nonlinearities  found  in  missile  autopilot 
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Figure  1:  Homing  Missile  Flight  Control  System 


design  problems,  including  saturation  limits,  coulomb  friction  and  backlash. 

Nonlinear  control  system  design  methods  will  be  described  in  Section  2.  These 
techniques  will  then  be  used  to  illustrate  autopilot  design  using  a  nonlinear  missile 
model  in  Section  3.  Conclusions  will  be  given  in  Section  4. 


2  Nonlinear  Autopilot  Design  Methods 

The  nonlinear  dynamic  models  used  for  autopilot  synthesis  are  assumed  to  be  of 
specified  in  the  form: 

x  =  f(x)+g  (x)  u  (1) 

Here,  x  is  the  state  vector,  u  is  the  control  vector,  f(x)  is  a  vector  of  state- 
dependent  nonlinear  functions  and  g{x)  is  a  matrix  of  state-dependent  nonlinear 
functions.  The  state  vector  components  for  autopilot  design  typically  consist  of 
angle  of  attack  a,  angle  of  sideslip  /3,  and  p,  q,  r,  the  pitch,  yaw  and  roll  body  rates. 
Note  that  the  nonlinear  functions  in  equation  (1)  could  be  in  the  form  of  lookup 
tables.  Every  nonlinear  control  design  technique  discussed  in  this  paper  assumes 
that  a  model  of  this  form  is  available.  This  model  is  termed  as  the  Design  Model 
in  the  following.  Although  this  model  may  include  nonlinearities,  it  should  not 
include  sensor  or  actuator  dynamics  in  order  to  ensure  that  their  internal  states 
are  not  used  in  autopilot  computations. 

In  most  missiles,  the  dynamics  may  contain  input  nonlinearities  such  as  satu¬ 
ration  and  cleadzone.  In  this  case,  the  model  may  actually  be  given  in  the  form: 

x  =  h(x,  u)  (2) 

Note  that  the  control  variables  appear  nonlinearly  in  expression  (2).  In  order 
to  facilitate  design  using  the  methods  described  in  this  paper,  the  model  in  (2) 
must  be  transformed  into  the  form  of  equation  (1).  This  can  be  accomplished  by 
introducing  a  new  set  of  control  variable  uc,  connected  to  the  system  through  a  set 
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of  input  dynamic  compensators.  The  input  compensators  can  be  of  any  desired 
form,  as  long  as  the  new  control  variables  appear  linearly  on  the  right  hand  side. 
For  instance,  input  compensators  can  be  in  the  form  of  pure  integrators,  such  that: 

x  =  h(x,  u),  u  =  uc  (3) 

Note  that  the  augmented  model(3)  is  now  in  the  standard  form  with  uc  as  the 
new  control  vector.  In  addition  to  allowing  transforming  the  missile  model  into 
the  standard  form,  the  designer  can  choose  the  dynamic  properties  of  the  input 
compensator  to  shape  the  frequency  content  of  the  signals  provided  to  the  actu¬ 
ators.  Selection  of  the  input  compensator  can  be  thought  of  as  another  degree  of 
freedom  available  to  the  in  the  autopilot  designer. 

Nonlinear  autopilot  design  techniques  discussed  in  this  paper  may  be  broadly 
classified  into  Transformation  Based  Methods  and  Direct  Methods.  This  clas¬ 
sification  is  based  on  the  way  the  nonlinear  design  techniques  utilize  the  system 
dynamic  model.  In  the  transformation-based  approaches,  the  given  dynamic  model 
is  first  transformed  either  to  the  Brunovsky  canonical  form  [1,  3,  21]  or  to  the  state- 
dependent  coefficient  form  [7-9].  Transformed  models  are  then  used  to  design  a 
feedback-linearized  autopilot  or  a  state-dependent  Riccati  equation  autopilot. 

Direct  Methods,  on  the  other  hand,  do  not  require  any  transformation  of  the 
given  nonlinear  system  model.  These  methods  employ  the  user-supplied  models 
in  the  standard  form  to  synthesize  the  controllers.  The  three  direct  design  tech¬ 
niques  discussed  in  this  paper  are:  a)  Quickest  Descent  method  [4],  b)  Recursive 
Back-Stepping  technique  [5]  and  c)  Predictive  Control  [6]  approach.  The  following 
subsections  will  describe  these  autopilot  design  techniques  in  further  detail. 

2.1  Feedback  Linearization  Method 

Feedback  linearization  techniques  have  been  used  for  flight  control  system  de¬ 
sign  for  over  two  decades.  This  technique  has  been  used  extensively  in  high- 
performance  aircraft  flight  control  system  design  [10  -  15]  and  for  missile  autopilot 
design  [16  -  19].  In  order  to  motivate  numerical  approach  to  the  feedback  lineariza¬ 
tion  method,  the  following  will  provide  a  brief  outline  of  the  technique.  Additional 
details  on  missile  autopilot  design  process  using  the  feedback  linearization  method¬ 
ology  can  be  found  in  References  16  and  17. 

The  first  step  in  the  feedback  linearization  approach  is  that  of  transforming 
the  system  dynamics  into  the  Brunovsky  [21]  canonical  form.  In  this  form,  the 
dynamic  system  under  consideration  is  in  the  form  of  decoupled  chains  of  inte¬ 
grators.  The  transformation  process  “pushes”  all  the  nonlinearities  in  the  system 
to  the  inputs,  thereby  enabling  the  construction  of  an  invertible,  state-dependent 
linearizing  map.  Next,  new  pseudo-control  variables  are  defined  as  the  product  of 
the  state-dependent  linearizing  map  and  the  actual  control  variables.  The  result¬ 
ing  dynamic  system  is  in  linear,  time  invariant  form  with  respect  to  pseudo-control 
variables.  If  the  pseudo-control  variables  are  known,  actual  control  variables  can 
be  computed  using  the  inverse  of  the  linearizing  map.  Since  the  state  variables  for 
computing  the  inverse  transform  are  obtained  from  feedback,  this  process  is  some¬ 
times  referred  to  as  global  linearization  using  feedback  or  feedback  linearization. 

Feedback  linearization  process  involves  the  selection  of  a  set  of  “leading  states” , 
which  are  then  repeatedly  differentiated  until  the  control  variables  appear  on  the 
right  hand  sides.  In  the  missile  autopilot  design  process,  the  angle  of  attack  a  and 
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the  angle  of  sideslip/3are  often  used  as  the  leading  states.  Differential  equations 
for  these  states  are  then  repeatedly  differentiated  until  the  fin  deflection  appears 
on  the  right  hand  sides.  For  instance,  the  fin  deflection  will  appear  in  the  second 
derivative  of  the  angle  of  attack: 

&=  fi{a,/3,P,q,r)  +  gu(a,/3)Sp  + g12(a,fi)6q  + g13(a,/3)5r  =  v3  (4) 

Here,  6p,  Sq,  Sr are  the  roll,  pitch,  and  yaw  fin  deflections.  State-dependent  non¬ 
linear  functions  /i,  <?n,  <712,513  generally  include  the  partial  derivatives  of  aerody¬ 
namic  forces  with  respect  to  angle  of  attack,  angle  of  sideslip  and  the  body  rates. 
At  any  given  values  of  the  state  vector,  these  functions  can  be  computed  numeri¬ 
cally  using  forward  difference  [22] . 

Direct  force  contributions  from  fin  deflections  are  generally  neglected  during 
this  transformation  process.  This  is  based  on  the  fact  that  the  lateral  acceler¬ 
ation  components  are  generated  primarily  by  changing  the  angle  of  attack  and 
angle  of  sideslip  using  pitch  rate  and  yaw  rate.  Pitch  and  yaw  rates  can  be  in¬ 
fluenced  through  the  moments  generated  by  the  fin  deflection.  Note  that  in  tail- 
controlled  missiles,  direct  forces  generated  by  fin  deflections  cause  the  well-known 
non-minimum  phase  response  of  the  autopilot. 

An  expression  for  angle  of  sideslip  can  also  be  similarly  derived. 

/?  =  f2(a,(3,p,q,r)  +  g2i{a,  fd)5p  +  g22{a,  (3)8q  +  g23{o:,P)Sr  =  v2  (5) 

The  feedback  linearization  process  is  analogous  to  the  transformation  of  linear 
dynamic  systems  into  the  controllable  canonical  form  [21].  Finally,  the  differential 
equation  for  roll  rate  may  be  given  as: 

P  =  h{&,  q ,  r)  +  g3i{a,  f3)5p  +  g32{a,  (3)5q  +  g33{a,  P)8r  =  v3  (6) 

The  variables  Vi,V2,v3  on  the  right  hand  sides  of  the  expressions  (4),  (5),  (6) 
are  the  pseudo-control  variables  in  the  pitch,  yaw  and  roll  axes.  Note  that  the 
system  is  linear  with  respect  to  the  pseudo-control  variables.  In  the  interests 
of  simplifying  notation,  Mach  number  and  altitude  dependencies  of  the  state- 
dependent  nonlinear  functions  on  the  right  hand  sides  of  the  expressions  (4)  -  (6) 
have  been  dropped.  However,  the  present  development  can  be  used  without  any 
change  to  include  that  case  also. 

After  the  system  is  transformed  into  the  feedback-linearized  form,  any  linear 
control  design  method  can  be  applied  to  derive  the  pseudo  control  variables.  This 
step  is  particularly  simple  because  the  feedback  linearized  dynamic  system  is  in 
the  form  of  decoupled  chains  of  integrators.  Techniques  such  as  pole  placement 
[23],  LQR  [24]  and  sliding  mode  control  [2]  can  be  directly  applied  to  derive  feed¬ 
back  control  laws  for  pseudo-control  variables.  For  instance,  using  either  the  pole 
placement  method  or  the  LQR  method,  the  pseudo-control  laws  for  regulating  the 
missile  dynamics  can  be  found  in  the  form: 


v\  =  I\\  a  +  K2  a 

(7) 

v2  =  K3p  +  KAp 

(8) 

II 

CO 

(9) 
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Once  the  pseudo-control  vector  is  available,  corresponding  fin  deflections  can  be 
computed  using  the  expressions  (4),  (5),  (6)  as: 


sp 

5n 

512 

513 

-1 

K±  a  +  K2  a  — 

h  ' 

Sq 

= 

521 

522 

523 

K3p+KAp- 

h 

Sr 

_  531 

532 

533 

_  Kbp-  f3 

The  foregoing  approach  can  be  modified  to  permit  lateral  acceleration  command 
tracking  by  augmenting  the  system  with  integral  tracking  error  states.  In  this 
case,  two  new  integral  error  states  are  defined,  with  the  corresponding  differential 
equations: 


hi(a,  (3, 5 pi  5q^  (5r) 

(11) 

^2(0^  /?,  5p,  5qi  (5r) 

(12) 

The  nonlinear  functions  hi  and  relate  the  angle  of  attack,  angle  of  sideslip  and 
the  fin  deflections  to  the  lateral  acceleration  components.  Expressions  (11)  and 
(12)  can  be  used  in  conjunction  with  the  angle  of  attack,  angle  of  sideslip  and 
body  rate  dynamics  to  derive  the  feedback  linearizing  transformations.  As  in  the 
case  of  regulator  design,  the  effect  of  fin  deflections  on  the  lateral  acceleration 
components  will  be  neglected  during  the  transformation  process. 

From  the  foregoing  discussions,  it  can  be  observed  that  the  numerical  imple¬ 
mentation  of  the  feedback  linearization  technique  requires  the  identification  of  the 
leading  states,  their  relationship  to  the  control  variables,  and  the  right-hand-sides 
of  the  system  dynamics.  This  information  can  be  used  by  a  numerical  differenc¬ 
ing  scheme  to  find  the  partial  derivatives  required  for  computing  the  feedback 
linearizing  transformations.  For  the  missile  autopilot  design  problem,  the  leading 
states  and  their  dominant  relationships  to  the  control  variables  can  be  symbolically 
expressed  as: 

Sp-^p  (13) 

5q  —>  q  ->  ol  — >  e*  (14) 

5r  — >  r  — >  (3  — >  ev  (15) 

Note  that  these  expressions  only  capture  the  dominant  relationships.  Coupling 
terms  in  the  transformation  can  be  computed  by  determining  the  dependence 
of  each  of  the  state  variables  on  other  state  variables.  Numbering  the  states 
ey,ez,a,  f3,p,q,r  sequentially,  and  the  control  variables  Sp,SqiSr  sequentially,  the 
symbolic  relationships  in  (13,  (14),  (15)  can  be  captured  in  the  form  of  a  matrix: 


5  0  0 

6  3  2 

7  4  1 


(16) 


Each  row  of  this  matrix  corresponds  to  a  control  variable,  and  each  column  corre¬ 
sponds  to  a  state.  For  instance,  the  first  row  suggests  that  the  roll  fin  deflection 
influences  the  roll  rate.  Similarly,  the  second  row  suggests  that  the  pitch  fin  de¬ 
flection  mainly  influences  the  pitch  rate,  which  in  turn  influences  the  angle  of 
attack.  The  angle  of  attack  strongly  influences  the  lateral  acceleration  component 
in  the  pitch  plane.  The  matrix  (16),  together  with  the  right-hand-sides  of  the 
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state  equations  can  used  to  configure  a  numerical  differencing  procedure  to  auto¬ 
matically  construct  feedback  linearizing  transforms  at  any  given  value  of  the  state 
vector. 

As  discussed  elsewhere  in  this  paper,  for  a  nonlinear  control  technique  to  be 
useful  in  applications,  it  should  be  able  to  directly  employ  a  simulation  model  of 
the  dynamic  system.  By  insisting  that  the  simulation  model  be  arranged  such 
that  the  integrators  are  included  in  a  separate  block  as  shown  in  Figure  2,  right- 
hand-sides  of  the  state  equations  can  be  readily  isolated  for  the  numerical  feedback 
linearization  process. 


I'sed  in 


Initial 

Conditions 


Figure  2:  Desired  Form  of  the  Simulation  Model 

A  software  package  is  currently  available  [25,  26],  that  uses  a  simulation  model 
of  the  form  shown  in  Figure  2  to  derive  feedback  linearizing  transforms,  and  the 
corresponding  inverse  transforms.  Numerical  algorithms  used  in  this  software  were 
developed  over  the  past  decade,  and  have  been  reported  in  References  13,  27  -  29. 
Once  the  model  is  feedback  linearized,  numerical  methods  for  linear  control  system 
design  [22]  can  be  employed  to  design  the  pseudo-control  laws. 

It  has  previously  been  observed  [11]  that  in  higher-order  dynamic  systems 
such  as  the  autopilot  design  problem,  the  nonlinear  controller  robustness  can  be 
significantly  enhanced  by  designing  multiple  time-scale  controllers.  Robustness 
in  time-scale  separated  controllers  result  from  the  fact  that  higher-order  partial 
derivatives  of  the  nonlinearties  on  the  right  hand  sides  of  the  nonlinear  model  are 
not  used  in  control  law  derivation.  Additionally,  time-scale  separated  nonlinear 
controllers  can  exploit  the  hierarchical  structure  of  the  system  states  to  simplify 
control  law  implementation.  Advantages  of  time  scale  separation  in  the  context 
of  missile  autopilot  design  have  been  previously  investigated  [17].  The  feedback 
linearization  procedure  outlined  in  the  foregoing  can  be  readily  adapted  to  allow 
for  time  scale  separation  of  the  system  dynamics.  In  this  case,  the  user  will  have 
the  additional  responsibility  for  identifying  the  state  variables  to  be  used  in  slow 
and  fast  time-scales. 

2.2  State-Dependent  Riccati  Equation  Method 

State  Dependent  Riccati  Equation  (SDRE)  method  [7  -  9]  is  another  technique 
that  uses  transformed  dynamic  model  for  nonlinear  controller  design.  By  defining 
the  nonlinear  dynamics  of  the  system  in  terms  of  state  dependent  matrices,  this 
technique  allows  the  derivation  of  nonlinear  controllers  using  techniques  similar 
to  that  of  the  LQR  [24]  technique.  The  first  step  in  the  SDRE  technique  is  the 
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transformation  of  the  user  specified  dynamic  model  given  in  the  standard  form 
x  =  f  (x)  +  g  (x)  u  into  the  State  Dependent  Coefficient  (SDC)  form  [7,  9]. 

x  =  A  (x)  x  +  g(x)  u  (17) 

The  matrix  A(x)  is  an  instantaneous  parameterization  of  the  state-dependent  non¬ 
linear  functions  f{x).  For  multivariable  systems,  infinite  number  of  such  realiza¬ 
tions  can  be  shown  to  exist  [7,  9].  However,  only  those  parameterizations  for  which 
the  pair  [A  ( x ) ,  g  (x)]  is  controllable  at  the  given  x  should  be  considered  for  the 
design. 

Note  that  the  SDC  parameterization  is  distinct  from  the  conventional  Taylor 
series  linearization.  Given  the  simulation  model  of  a  dynamic  system  in  the  form 
illustrated  in  Figure  2,  instantaneous  SDC  parameterization  can  be  obtained  by 
evaluating  the  vector  nonlinear  function  f(x)  using  a  set  of  linearly  independent 
probe  vectors  (2.  . ,  C «•  As  a  practical  matter,  since  the  behavior  of  the  nonlin¬ 

earities  in  the  neighborhood  of  the  current  system  state  are  not  explicitly  known, 
it  is  wise  to  choose  probe  vectors  that  are  close  to  the  current  state  vector.  At  the 
current  state  vector,  the  vector  nonlinear  function  /( x)  can  be  extracted  from  the 
simulation  model  by  setting  the  control  vector  to  zero. 

The  probe  vectors  are  constructed  by  adding  small  magnitude  perturbation 
vectors  crysg(Tzsg<74Suuuuc7„sgoghe  nominal  state  vector  to  yield  a  set  of  linearly 
independent  vectors: 

C2  =  X  +  <T2,  <2  =  X  +  <73 ,  C3  =  X  +  CT4,  . ,  Cra  =  X  +  On  (18) 

The  nonlinear  function  f{x)  is  next  evaluated  using  these  linearly  independent 
vectors  to  assemble  a  matrix  equation  of  the  form: 

[  f(x)  /(C 2)  .  /( Cn)  ]  =  A  (x)  [  X  C 2  .  C n  ]  (19) 

At  any  given  value  of  x,  this  linear  matrix  equation  can  be  solved  for  the  elements 
of  A(x).  Since  the  probe  vectors  and  the  state  vector  are  linearly  independent, 
this  equation  is  well  conditioned,  and  can  be  solved  using  well-known  numerical 
linear  algebraic  algorithms. 

Note  that  the  foregoing  computations  will  have  to  be  carried  out  at  every 
sample.  The  SDC  matrix  A(x)  from  these  computations  can  next  be  used  to 
formulate  and  solve  the  SDRE  control  problem. 

As  an  aside,  it  is  interesting  to  examine  the  relationship  between  the  numerical 
construction  of  the  SDC  model  and  the  conventional  Taylor  series  linearization.  If 
the  perturbation  vectors  crysgcrzsgcr4suuuucr„sgre  small,  it  can  be  found  that: 

df 

— ,  at  x  =  0  (20) 

ox 

Note  that  this  corresponds  to  the  Taylor  series  linearization  of  the  system  dynam¬ 
ics  about  the  origin.  Thus,  the  present  numerical  methodology  for  constructing 
the  SDC  model  automatically  reverts  to  Taylor  series  linearization  of  the  system 
dynamics  near  the  origin  of  the  state  space.  For  constant  control  influence  ma¬ 
trix  case,  the  present  SDC  parameterization  scheme  preserves  the  controllability 
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properties  of  the  dynamic  system  near  the  origin.  Since  the  only  restriction  on  the 
probe  vectors  is  that  they  be  linearly  independent,  it  is  possible  to  construct  an 
infinite  variety  of  SDC  parameterizations  for  a  given  dynamic  system. 

The  SDC  form  of  the  system  state  equations  is  next  used  to  cast  the  control 
problem  as  an  infinite-horizon  nonlinear  regulator  minimizing  the  cost: 


J  = 


J  xtQ(x)x  +  uT*(*)u  dt(21) 


*0 


subject  to  the  nonlinear  differential  constraint:  x  =  A  ( x )  x  +  g  ( x )  u.  The  state- 
dependent  matrices  Q(x)  and  R(x )  are  chosen  by  the  designer  to  achieve  desired 
properties  of  the  closed-loop  system. 

References  7  and  9  show  that  the  solution  to  this  problem  can  be  obtained  by 
solving  a  state  dependent  Riccati  equation: 

AT  (x)  P(x)  +  P  (x)  A(x)  —  P  (x)  g  (x)  R~x  (x)  gT  (x)  P  (x)  +  Q  =  0  (22) 


The  state  dependent  feedback  gain  can  then  be  computed  as: 

K  (x)  =  i?_1  (x)  gT  (x)  P  (x)  (23) 

The  SDRE  nonlinear  control  law  is  of  the  form: 

u  =  — K  (x)  x  (24) 


Reference  7  has  shown  that  under  rather  mild  restrictions  on  Q(x)  and  R(x),  the 
SDRE  control  law  will  globally  stabilize  the  nonlinear  dynamic  system.  Figure 
3  illustrates  the  computational  steps  involved  in  the  SDRE  technique.  At  each 
time  step,  SDC  parameterization  of  the  dynamic  system  is  generated  and  used 
to  formulate  the  algebraic  Riccati  equation.  This  equation  is  then  solved  for  the 
state-dependent  feedback  gains. 


Off-Line  Selection  of 
State  &  Control  | 
Weighting  Matrices 
Q(x),  R(x)  4 


State 
Vector:  x 
From  the 
Dynamic  Systeml 


Compute 
Model  Matrices 

A(x),  B(x) 

Formulate  and  Solve 
The  StateDependent 
Riccati  Equation  for  P(x)  >0 


Compute 

State-Dependent  Gain 
K  (x)  =  R-l(x)  Bt(x)  P(x) 

y 


Control 
Vector:  u 
To  the 


State  Feedback 
•  Control  Law: 
u  =  -  K(x)  x 


Dynamic  System 


Figure  3:  SDRE  Control  Computations 


Note  that  the  SDRE  technique  is  computationally  demanding,  requiring  the 
solution  of  a  7x7  algebraic  Riccati  equation  at  each  sample.  However,  previous 


182 


P.K.  Menon  ,  E.  J.  Ohlmeyer 


research  [30]  has  demonstrated  that  these  computations  are  well  within  the  capa¬ 
bilities  of  commercially  available  processors. 

2.3  Quickest  Descent  Method 

The  Quickest  Descent  method  discussed  in  Reference  4  is  a  Lyapunov  function 
optimizing  approach  to  nonlinear  feedback  controller  design.  This  is  perhaps  the 
simplest  of  all  the  five  nonlinear  control  techniques  discussed  in  this  paper.  In  this 
approach,  the  control  problem  is  viewed  as  a  function  minimization  problem  in 
the  state  space.  A  descent  function,  W(x)  satisfying  certain  specified  properties 
is  first  selected.  Control  vector  is  then  chosen  to  minimize  this  descent  function. 

The  descent  function  W(x)  is  required  to  be  bounded,  continuous  and  con¬ 
tinuously  differentiable  in  the  region  of  interest.  In  addition,  the  target  state  is 
required  to  be  contained  within  the  region  of  interest.  Note  that  these  require¬ 
ments  are  more  restrictive  than  the  choice  of  Lyapunov  functions.  Although  no 
general  guidelines  are  available  for  the  selection  of  descent  functions,  it  appears 
that  physical  quantities  such  as  the  total  energy  in  the  system  can  be  used  as 
the  starting  point.  In  problems  such  missile  autopilot  design,  a  properly  defined 
quadratic  function  of  the  states  can  be  used  as  the  descent  function. 

Once  the  descent  function  W(x)  is  selected,  the  feedback  control  u(x)  is  chosen 
so  that  W(x)  decreases  at  each  state  of  the  system.  If  the  minimization  process 
is  cast  as  a  steepest  decent  optimization  problem,  the  resulting  technique  can  be 
termed  as  the  steepest  decent  control  methodology.  Reference  4  shows  that  a 
more  direct  approach  is  to  choose  the  control  variables  to  minimize  the  time-rate 
of  change  of  the  descent  function  approach.  In  this  case,  the  control  methodology 
can  be  termed  as  the  Quickest  Decent  method. 

In  the  quickest  descent  method,  control  is  obtained  by  minimizing  the  time 
rate  of  change  of  the  descent  function.  Thus,  the  optimization  problem  is  of  the 
form: 


min 

U 


dW  ( x ) 
dt 


or 


min 

u 


'dW 

dx 


{/  0)  +  9  0)  u} 


(25) 


Since  the  control  variable  appears  linearly  in  the  system  dynamics,  the  optimiza¬ 
tion  problem  is  meaningful  only  if  the  control  variables  are  constrained.  The 
control  constraints  can  be  specified  in  the  form:  |u|  <  wmax.  The  quickest  descent 
control  is  then  given  as: 

f  dW  ) 

If  <  -Q^g  (x)  J-  >  0,  u  =  umin  (26) 

f  dW  ) 

If  \~g^r9{x)>  <  o,  u  =  umax  (27) 


Note  that  the  control  is  bang-bang.  Under  the  present  formulation  of  the  quick¬ 
est  descent  method,  the  control  variables  will  chatter  between  their  limits  as  the 
system  approaches  the  minimum  of  the  descent  function  W (a:) . 

In  missile  autopilot  design  problem,  the  descent  function  can  often  be  specified 
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as  a  quadratic  function  of  the  form: 

e* 
ey 
a 

P 

P 

q 

r 

with  P  being  a  positive  definite  matrix.  The  elements  of  the  matrix  P  has  to  be 
carefully  to  ensure  that  the  control  objectives  are  satisfied.  For  instance,  diagonal 
elements  of  P  can  be  chosen  to  drive  state  vector  components  such  as  the  integral 
tracking  error  states  and  the  body  rates  to  zero,  while  the  off-diagonal  elements  can 
be  chosen  to  preserve  the  coupling  between  the  missile  state  variables  to  achieve 
the  desired  response.  Selection  of  the  control  bounds  and  the  matrix  P  are  the 
two  degrees  of  freedom  available  to  the  designer  in  the  quickest  descent  control 
technique. 

2.4  Recursive  Back-stepping  Method 

As  can  be  observed  from  Section  2.1,  the  feedback  linearization  method  can¬ 
cels  the  system  nonlinearities  and  replaces  them  with  a  linear  dynamic  system. 
The  main  premise  behind  the  recursive  back  stepping  technique  [5]  is  that  certain 
portions  of  the  system  nonlinearities  are  worth  preserving.  This  objective  is  sat¬ 
isfied  by  formulating  the  control  problem  using  the  second  method  of  Lyapunov. 
However,  since  there  is  no  direct  approach  for  constructing  Lyapunov  functions  for 
multivariable  nonlinear  dynamic  systems,  the  backstepping  procedure  relies  on  a 
recursive  procedure.  Just  as  the  system  nonlinearities  were  pushed-back  to  the  in¬ 
puts  in  the  feedback  linearization  methodology,  recursive  back  stepping  technique 
constructs  Lyapunov  function  for  the  nonlinear  dynamic  system  by  stepping-back 
from  the  output  state  variables  to  the  controls.  In  some  respects,  this  technique 
bears  a  strong  resemblance  to  the  multiple  time-scale  feedback  linearization  design 
technique. 

The  recursive  back  stepping  design  technique  assumes  that  the  model  is  spec¬ 
ified  in  a  triangular  form  as  shown  in  equations  (29)  -  (31): 

xi  =  fi  (xi)  +  gi  (xi)  x2  (29) 

X2  =  fi  (xi,x2)  +  gi  (xi,x2)  x3  (30) 

Xn  =  fn{xi,x2,---  ,xn)  +  gn{xi,x2,  -  ■■  ,x„)  u  (31) 

Here,  x\,x2,.  ■  ■  xn&re  the  components  of  the  state  vector.  Each  scalar  system  is 
stabilized  with  the  following  state  as  the  control  variable.  For  example,  observes 
as  the  control  variable  for  X\-  dynamics,  x%  for  x2-  dynamics  and  so  on.  Note  that 
such  a  structure  can  be  found  in  missile  dynamics,  if  the  direct  force  contributions 
arising  from  fin  deflections  are  neglected.  The  triangular  structure  in  missile  dy¬ 
namics  consists  of  the  pitch  and  yaw  rates  generating  angle  of  attack  and  angle  of 
sideslip,  which  in  turn  result  in  the  lateral  acceleration  components. 


W  (x)  =  ez  ey  a  (3  p  q  r  P 
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Controllers  are  synthesized  for  each  scalar  dynamic  system  using  the  second 
method  of  Lyapunov.  At  each  stage  of  the  backstepping  process,  Lyapunov  func¬ 
tions  are  selected  to  preserve  desired  nonlinearities.  In  most  practical  problems, 
quadratic  Lyapunov  functions  can  be  used  to  realize  most  of  the  benefits  of  this 
method. 

In  the  missile  autopilot  example,  the  back  stepping  process  begins  with  the 
integral  tracking  error  states  ez,  ey  and  the  roll  rate  p.  Control  system  design  for 
the  roll  channel  does  not  involve  any  backstepping,  since  the  roll  fin  deflection  is 
directly  related  to  the  roll  acceleration.  A  quadratic  Lyapunov  function  of  the 
form:  |p2 can  be  used  to  derive  the  roll  fin  deflection  as  a  function  of  the  missile 
states  that  will  drive  the  roll  rate  to  zero. 

In  the  pitch  and  yaw  channels,  angle  of  attack  and  angle  of  sideslip  form  the 
control  variables  for  the  integral  tracking  error  states  in  the  first  stage  of  the  back 
stepping  process.  A  Lyapunov  function  of  the  form:  |  [e2  +  e2]  can  be  used  to 
derive  the  angle  of  attack  a±  and  angle  of  sideslip  /3i  that  will  drive  the  integral 
tracking  error  states  to  zero.  Next,  the  Lyapunov  functions  are  augmented  by 
quadratic  terms  in  (a.\  -ap and  (/3i~/3p  to  form  the  second  stage  of  the  backstepping 
process.  In  the  second  stage,  the  pitch  rate  q2  and  the  yaw  rate  r2  are  the  “control 
like”  variables.  These  variables  are  chosen  to  ensure  that  the  time  rate  of  change 
of  the  augmented  Lyapunov  function  will  be  less  than  or  equal  to  zero.  In  the 
last  stage  of  the  back  stepping  process,  the  Lyapunov  function  is  augmented  by 
quadratic  terms  in  ((72  —  q )  and  (r^  —  r).  Pitch  and  yaw  fin  deflections  that  drive 
the  time  rate  of  change  of  the  resulting  Lyapunov  function  to  be  less  than  or  equal 
to  zero.  At  each  stage  of  the  backstepping  process,  the  intermediate  states  are 
eliminated  using  the  given  state  equations,  see  Reference  5  for  details. 

In  numerical  implementations  of  the  backstepping  process,  the  user  can  specify 
the  backstepping  sequence  through  a  matrix  similar  to  that  used  for  specifying  the 
feedback  linearization,  expression  (15).  The  user  can  also  specify  the  form  of  the 
Lyapunov  functions  to  be  used  in  each  step  of  the  backstepping  procedure.  Since 
the  nonlinearities  contained  in  the  missile  dynamics  is  generally  limited  to  lookup 
tables,  products  of  the  state  variables  and  transcendental  functions,  quadratic 
Lyapunov  functions  are  often  adequate  to  obtain  good  response  from  the  autopilot. 
This  information  together  with  a  simulation  model  of  the  form  shown  in  Figure  2 
are  sufficient  to  develop  a  numerical  procedure  for  recursive  backstepping  process. 

2.5  Predictive  Control 

The  predictive  control  methodology  [6]  has  been  popular  in  the  process  control 
industry  for  the  past  two  decades.  In  this  technique,  a  control  history  that  will 
drive  the  system  states  to  the  desired  values  at  the  end  of  a  prediction  interval 
is  computed  at  every  sample  interval  using  an  optimization  algorithm.  Most  of 
the  techniques  described  in  the  literature  employ  multi-step  predictors  to  imple¬ 
ment  the  controllers.  In  nonlinear  systems,  the  use  of  predictive  control  technique 
requires  the  use  of  an  on-line  iterative  optimization  algorithm.  Since  the  conver¬ 
gence  of  optimization  techniques  for  general  nonlinear  problems  is  not  assured, 
the  performance  of  the  predictive  control  technique  cannot  be  guaranteed  in  these 
situations 

However,  if  the  nonlinear  dynamic  system  is  of  the  standard  form  given  by 
equation  (1),  and  if  the  control  problem  is  cast  as  a  one-step-ahead  predictive 
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control,  the  system  performance  becomes  more  predictable.  This  approach  is 
adopted  in  the  present  formulation  of  the  nonlinear  predictive  control  technique. 
Thus,  the  control  problem  is  cast  as  the  minimization  of  the  quadratic  objective 
function 

JJ±x[+1  Q(x,t)  xi+1  +  uf  R  (x,t)  11^32) 

with  respect  to  the  control  variable  u,  subject  to  the  differential  constraint: 

x  =  f(x)+g  (x)  u  (33) 

The  differential  constraint  can  be  used  to  eliminatea:,;_|_i  from  the  objective  function 
by  defining  a  numerical  integration  algorithm.  Any  one  of  the  several  numerical 
integration  techniques  can  be  used  for  this  purpose.  However,  if  the  nonlinear 
dynamic  system  is  specified  in  the  standard  form,  the  linearity  of  the  differential 
constraints  with  respect  to  the  control  variables  can  be  preserved  if  techniques 
such  as  Euler’s  integration  method  or  the  Adams-Bashforth  numerical  integration 
scheme  [22]  are  employed.  For  instance,  if  the  Euler’s  integration  formula  is  used 
to  integrate  the  differential  constraint  (33)  with  a  step  size  At,  the  control  vector 
that  minimizes  (32)  can  be  found  to  be: 

uk  =  -  [R-1  +  At2  g[  Q  gk\  1  [R~x  +  At  gl  Qxk  +  At2g%Q  fk]  (34) 

The  use  of  multi-step  integration  formulae  will  produce  more  complex  expressions 
for  one-step-ahead  optimal  control.  The  state  and  control  weighting  matrices  Q 
and  R  must  be  chosen  to  satisfy  the  descent  property: 

Ji+i  <  J;  (35) 

As  in  the  Quickest  Descent  method,  the  structure  of  the  state  weighting  matrix 
must  be  chosen  to  establish  the  desired  relationships  between  the  state  variables. 


3  Design  Example  —  Nonlinear  Autopilot  Design 
Using  Missile  Longitudinal  Dynamics 

The  application  of  the  nonlinear  control  system  design  methods  described  in  Sec¬ 
tion  2  will  be  applied  to  the  missile  autopilot  design  problem  in  this  section.  In  the 
interests  of  a  compact  presentation,  formulation  of  the  nonlinear  regulation  prob¬ 
lem  will  only  be  examined.  As  discussed  in  Section  2,  extension  to  the  command 
tracking  case  involves  the  augmentation  of  the  system  dynamics  with  integral 
tracking  error  states,  and  is  direct.  Longitudinal  dynamic  model  of  a  missile  used 
to  illustrate  the  nonlinear  autopilot  design  techniques  in  this  paper  is  obtained 
from  Reference  8. 

The  dynamics  of  a  generic  homing  missile  coasting  in  the  vertical  plane  is  given 
by  the  expressions  (36)  through  (37).  The  model  incorporates  Mach  number  M, 
angle  of  attack  a,  flight  path  angle  7  and  pitch  rate  q  as  the  state  variables  and 
the  pitch  fin  deflection  S  as  the  control  variable.  The  aerodynamic  axial  force, 
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normal  force  and  the  pitching  moment  have  all  been  expressed  in  terms  of  the 
state  variables  on  the  right  hand  sides  of  these  differential  equations. 


M  =  0.4008  M2o?  sin  (a)  -  0.6419  M2|a|Qsin(a)"0-2010 

9  n  rum  sin(“)  a_0-0311  sin(7) 

-0.0062  m2-0'0403  M 


M2(2-T)“sin(“> 


(36) 


a  =  0.4008  Mo?  cos  (a)  —  0.6419  M  |a|  a  cos  (a)  —  0.2010  M  (2  —  4^)  a  cos  (a) 

-0.0403  Mcos  (a)  (5+0.0311  +  q 

(37) 


7  =  —0.4008  Mo?  cos  (a)  +  0.6419  M  |a|  a  cos  (a)  +  0.2010  M  (2 
+0.0403  Mcos  (a)  5-0.0311 


=y)  a  cos  (a) 


(38) 


q  =  49.82  M2a3  -  78.86M 


2|a|a-f-3.60  M2(  7 


M2 


(39) 


-2.12M29 


Missile  autopilot  design  problem  is  concerned  with  the  short-period  dynamics  with 
a  and  q  as  the  state  variables.  The  flight  path  angle  7  and  Mach  number  M,  are 
necessary  for  calculating  the  aerodynamic  and  gravitational  forces.  The  differential 
equations  for  7  and  M  can  be  treated  as  auxiliary  expressions  that  contribute  to  the 
right  hand  sides  of  the  short  period  dynamics.  Pitch  fin  deflection,  <5  is  the  control 
variable.  Note  that  the  gravitational  acceleration  term  is  normally  neglected  in 
missile  autopilot  design  problems. 

As  the  first  step  in  the  autopilot  design  process,  a  computer  simulation  of  the 
missile  is  constructed  using  the  system  dynamics.  A  step  size  of  10ms  was  used 
in  the  simulations.  The  pitch  rate  is  identified  as  the  first  state,  and  the  angle 
of  attack  is  denoted  as  the  second  state  in  the  system.  Both  these  states,  Mach 
number  and  the  flight  path  angle  were  assumed  to  be  available  from  measurements. 
No  actuator  dynamics  was  included  in  the  simulation.  Design  parameters  and 
simulation  results  for  each  of  the  five  techniques  will  be  given  in  the  following 
sections. 

3.1  Nonlinear  Regulator  Designs  Using  Feedback  Linearization 

The  first  step  in  this  design  technique  is  that  of  generating  the  feedback  lin¬ 
earization  map.  In  the  autopilot  design  example,  the  feedback  linearization  se¬ 
quence  can  be  specified  as: 

<5  — >  q  — >  a  (40) 

Since  the  pitch  rate  and  angle  of  attack  states  have  been  sequentially  numbered, 
this  can  also  be  represented  by  a  row  vector:  [1  2].  This  notation  implies  that  the 
fin  deflection,  the  only  control  variable  in  the  problem,  will  be  used  to  generate 
pitch  rate,  which  would  then  produce  the  desired  angle  of  attack.  Thus,  the 
feedback  linearization  map  will  be  constructed  by  first  differentiating  the  right 
hand  side  of  the  differential  equation  for  angle  of  attack,  followed  by  a  substitution 
of  the  right  hand  side  of  the  pitch  rate  equation.  Forward  difference  method 
was  used  for  the  numerical  computation  of  the  partial  derivatives.  The  state 
perturbations  used  for  these  computations  were  10~4Radians. 
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Feedback  linearized  form  of  the  system  dynamics  is: 


(41) 


This  model  can  be  used  to  design  feedback  control  laws  for  the  pseudo-control 
variable.  The  control  law  can  then  be  inverse  transformed  to  obtain  fin  deflections 
that  will  regulate  the  nonlinear  dynamic  system.  This  process  will  be  illustrated 
using  three  different  control  approaches. 

(a)  Pole  Placement  design: 

The  poles  of  the  feedback  linearized  dynamic  system  were  chosen  to  provide 
critical  damping  and  a  natural  frequency  of  about  14  rad/sec.  The  response  of  the 
nonlinear  missile  dynamics  under  this  feedback  linearized  control  law  is  illustrated 
in  Figures  4  and  5.  It  can  be  observed  that  the  system  responses  are  essentially 
that  of  the  feedback  linearized  system.  The  responses  have  been  found  to  re¬ 
main  invariant  even  under  ±10%  perturbations  in  aerodynamic  force  and  moment 
models. 

(b)  Linear  Quadratic  Regulator  design: 

The  LQR  approach  is  used  next  to  design  the  feedback  linearized  control  law. 
The  state-weighting  matrix  was  chosen  to  be: 


1000  0 
0  1 


(42) 


The  control  weighting  was  chosen  as  0.01. 

The  angle  of  attack  and  pitch  rate  responses  to  the  initial  conditions  are  illus¬ 
trated  in  Figure  6.  Corresponding  fin  deflection  time  history  is  given  in  Figure  7. 
It  may  be  observed  that  the  responses  from  the  feedback  linearized  LQR  design  are 
comparable  to  the  ones  from  the  pole  placement  design.  As  with  that  case,  ±10% 
perturbations  in  the  system  model  had  no  observable  effects  on  the  closed-loop 
system  responses. 

(c)  Sliding  Mode  design: 

Sliding  mode  control  methodology  from  Reference  2  is  next  used  to  design  a 
regulator  for  the  feedback  linearized  system  dynamics.  The  sliding  surface  param¬ 
eter  A  is  chosen  as - 20,  uncertainty  parameters  are  set  to  zero,  the  convergence 

parameter  to  the  sliding  surface  ?y  is  chosen  as  20,  and  the  thickness  of  the  bound¬ 
ary  layer  e  is  chosen  as  0.1.  The  time  histories  of  the  angle  of  attack  and  pitch 
rate  for  the  feedback  linearized  sliding-mode  autopilot  are  given  in  Figure  8.  Cor¬ 
responding  fin  deflection  history  is  in  Figure  9.  Note  that  the  response  of  the 
closed-loop  system  is  significantly  different  from  that  of  the  pole  placement  or 
LQR  designs. 

3.2  State  Dependent  Riccati  Equation  Autopilot  Design 

State  perturbations  used  for  constructing  the  SDC  form  of  the  missile  model 
were  10_4Radians.  Although  the  design  technique  allows  state-dependent  state 
and  control  weighting  matrices,  constant  weights  were  used  in  the  present  case. 
These  were: 


Q 


100  0 
0  10  ’ 


R  =  10(43) 
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Figure  4:  Angle  of  Attack  and  Pitch  Rate  Responses  for  the  Feedback  Linearized 
Regulator  with  Pole  Placement  Design 
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Figure  5:  Fin  Deflection  for  the  Feedback  Linearized  Regulator  with  Pole  Place¬ 
ment  Design 


Figure  6:  Angle  of  Attack  and  Pitch  Rate  Responses  for  the  Feedback  Linearized 
Regulator  with  LQR  Design 
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Figure  7:  Fin  Deflection  for  the  Feedback  Linearized  Regulator  with  LQR  Design 


Angle  of  Attack  vs.  Time 


Pitch  Rate  vs.  Time 


Time  (sec) 


Figure  8:  Angle  of  Attack  and  Pitch  Rate  Responses  for  the  Feedback  Linearized 
Regulator  with  Sliding  Mode  Design 
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Figure  9:  Fin  Deflection  for  the  Feedback  Linearized  Regulator  with  Sliding  Mode 
Design 


Response  of  the  SDRE  autopilot  are  given  in  Figures  10  and  11.  Note  that  it 
is  possible  to  further  improve  the  speed  of  response  by  choosing  alternate  design 
weights. 

3.3  Quickest  Descent  Design 

The  design  parameters  for  the  Quickest  Descent  technique  are  the  descent 
function  and  the  limits  on  the  control  variables.  As  discussed  elsewhere  in  this 
paper,  the  descent  function  can  be  chosen  as  a  quadratic  function  in  most  practical 
design  problems.  For  the  autopilot  design,  the  descent  function  was  chosen  to  be 
of  the  form: 


W  (x)  =  a  q  P 


(43) 


The  weighting  matrix  and  the  control  limits  were  chosen  to  be: 


5000  0.5 
0.5  1 


\u\  <  0.01 


(44) 


The  state  and  control  histories  obtained  from  a  closed-loop  simulation  of  the  quick¬ 
est  descent  autopilot  are  given  in  Figures  12  and  13.  It  may  be  observed  that  the 
control  chatters  after  an  initial  bang-bang  region. 

3.4  Recursive  Backstepping  Autopilot  Design 

Recursive  backstepping  autopilot  is  synthesized  using  the  angle  of  attack  and 
pitch  rate  as  the  backstepping  variables,  and  fin  deflection  as  the  control  variables. 
Numerical  partial  derivatives  required  in  the  backstepping  procedure  were  com¬ 
puted  using  forward  difference  technique  [2]  with  10_4Radians  perturbation  in  the 
states. 

Quadratic  Lyapunov  functions  were  used  in  each  step  of  the  process,  with  a 
convergence  rate  of  20  in  each  stage.  The  angle  of  attack  and  pitch  rate  time 
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Figure  10:  Angle  of  Attack  and  Pitch  Rate  Responses  for  the  Feedback  Linearized 
Regulator  with  SDRE  Design 


Figure  11:  Fin  Deflection  for  the  Feedback  Linearized  Regulator  with  SDRE  De¬ 
sign 
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Pitch  Rate  vs.  Time 


Figure  12:  Angle  of  Attack  and  Pitch  Rate  Responses  for  the  Feedback  Linearized 
Regulator  with  Quickest  Descent  Design 


Figure  13:  Fin  Deflection  for  the  Feedback  Linearized  Regulator  with  Quickest 
Descent  Design 
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histories  for  the  recursive  backstepping  autopilot  are  shown  in  Figure  14.  Fin 
deflection  history  is  given  in  Figure  15. 


Angle  of  Attack  vs.  Time 


Pitch  Rate  vs.  Time 


Figure  14:  Angle  of  Attack  and  Pitch  Rate  Responses  for  the  Backstepping  Au¬ 
topilot  Design 


3.5  Predictive  Autopilot  Design 

Predictive  autopilot  design  used  a  second-order  Adams-Bashforth  integration 
algorithm  to  formulate  the  predictive  performance  index.  The  state  weighting 
matrix  and  the  control  weights  were  chosen  as: 


Q 


10  0 
0  10  ’ 


R=  2 


Note  that  this  design  equally  weights  angle  of  attack  response  and  the  pitch 
rate  response.  The  state  vector  and  control  response  of  the  predictive  autopilot  are 
illustrated  in  Figures  16  and  17.  As  with  other  design  techniques  with  quadratic 
criteria,  the  speed  of  response  of  the  predictive  autopilot  can  be  improved  by 
reducing  the  control  weighting  or  by  increasing  the  state  weights. 


4  Conclusions 

This  paper  discussed  five  different  approaches  for  computer-aided  nonlinear  autopi¬ 
lot  design.  These  are  numerical  implementations  of  the  nonlinear  control  system 
design  techniques  discussed  in  the  literature.  Consequently,  they  can  use  numerical 
simulation  models  of  the  dynamic  systems  of  arbitrary  complexity  to  automati¬ 
cally  derive  stable  nonlinear  closed-loop  control  systems.  The  focus  of  the  present 
paper  was  on  the  application  of  these  techniques  to  the  design  missile  autopilots. 
General  design  approaches  were  first  outlined,  and  illustrated  using  a  nonlinear 
model  of  a  missile.  The  approaches  presented  in  this  paper  have  been  employed  to 
design  full-order  autopilots  and  integrated  guidance-control  systems  using  a  vari¬ 
ety  of  missile  models.  The  examples  presented  in  this  paper  show  that  numerical 
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Figure  15:  Fin  Deflection  history  for  the  Backstepping  Autopilot  Design 


Angle  of  Attack  vs.  Time 


Pitch  Rate  vs.  Time 
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Figure  16:  Angle  of  Attack  and  Pitch  Rate  Responses  for  the  Predictive  Autopilot 
Design 
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Figure  17:  Time  history  of  the  Fin  Deflection  for  the  Predictive  Autopilot  Design 


approach  to  nonlinear  control  system  design  is  feasible,  and  can  be  carried  out 
with  a  level  of  confidence  comparable  to  that  of  linear  design  techniques. 
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